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Abstract 

The isotropic rigid and non-rigid rotators in the framework of Tsallis statistics 
are studied in the high and low temperature limits. The generalized partition 
functions, internal energies and heat capacities are calculated. Classical re- 
sults of the Boltzmann-Gibbs statistics have been recovered as non-extensivity 
parameter approaches to 1. It has also been observed that non-extensivity 
parameter q behaves like a scale parameter in the low temperature regime of 

the rigid rotator model. 
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I. INTRODUCTION 



The non-extensive generalization of the standard Boltzmann-Gibbs statistics was pro- 
posed in 1988 by C. Tsallis[l-4]. This non-extensive generalization begins with the supposi- 
tion of a new, fractal inspired entropy S q = k{l — J2iPi)/(q ~ 1) where q is any real number 
and k is a positive constant which becomes Boltzmann constant in the limit 1. In this 
case, the new entropy S q takes the form -ke J2iPi In Pi, which is simply the standard entropy 
formula. This new non-extensive statistics has been studied a great deal to list its properties 
and it is also applied to many well-known examples of the Boltzmann-Gibbs statistics such 
as Ehrenfest theorem[5], Fokker- Planck equations[6] and quantum statistics [7]. 

We first calculate the partition function of the isotropic rigid rotator for both high 
energy and low energy limits and then obtain the generalized internal energy and specific 
heat expressions for both cases. This forms the Section II. The usual Maxwell-Boltzmann 
results are obtained in the limit q— > 1. In Section III, we turn our attention to non-rigid 
rotator and follow the same procedure for it as in the Section II. Results and discussions are 
given in Section IV. 



where j=0,l,2, /j, = m\mij{m\ + m 2 ) reduced mass of the nuclei and a being the equilib- 
rium distance between them. The degeneracy of each level is gj — 2j + 1. In non-extensive 
formalism (NEXT), the partition function is given by 



II. ISOTROPIC RIGID ROTATOR 



Energy levels of a rigid rotator are [8] 



E 3 = 3 ( 3 + l)h 2 /2^a 2 . 
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Now we are ready to look for its analytic solutions in the high and low temperature limits. 
i)High temperature limit 

At high temperatures, j, <C 1 and [l_ilz^Mi±lM] 1/(1-9) changes slowly as j changes. So 
we take it as a continuous function of j. Letting j(j+l)=x and substituting (2j+l)dj=dx, 
then, we have 
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For the interval 1 ( q ( 2, solution to the above integral is 



As q goes to 1, we get 
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Z,-i = j- (6) 



The expression above is exactly the Maxwell-Boltzmann (MB) partition function for the 
isotropic rigid rotator in the high temperature limit. Thus, we calculate the generalized 
internal energy function from 

Uq = « 32zl, (7) 

dp 1 — q 

in the non-extensive case. It becomes 

where a = A. It is important to see that again, in the limit when q approaches to 1, we 
obtain the result already known in MB statistics, i.e., 

£W = ^ = k B T. (9) 

It is also possible to calculate the specific heat in this statistics as 



c, - f£, (10) 
where XJ q is the internal energy function. We immediately get 

^ = (^) 1 -^( 2 -^- (11) 

As q — ► 1, C q — > k B is verified easily. All these calculations are carried out by using the 
constraint Y^LiPi £ i — U q . But, this choice of internal energy constraint presents several 
difficulties. For example, the probability distribution obtained by using this constraint is not 
invariant under uniform translation of the energy spectrum. Another important consequence 
of this constraint is that it leads to violation of energy conservation macroscopically.The 
solution to all these difficulties are recently proposed by C. Tsallis et al. [9]. They introduced 
the new internal energy constraint as follows 

= U q . (12) 

This choice of constraint remedied all previous difficulties. Since we have studied all ther- 
modynamical quantities with the previous internal energy constraint, we must reconsider 
them with this new constraint. There are two ways to do this: Firstly, we can recalculate 
all thermodynamical quantities with this new internal energy constraint by forming the new 
partition function. Another method is to find the relation between the temperatuire param- 
eters of the old and new calculations. If one has all thermodynamical functions calculated 
with the old constraint and the relation between temperature parameters is known, it is 
possible to modify all previous calculations carried out with the constraint Y^LiPi £ i — U q . 
This is the method we will follow, because we already have the solutions with respect to old 
constraint. 

We begin by writing all previous calculations in terms of intermediate variable t' where 
t' = l/(/3'e). From now on, The superscript (2) will refer to calculations done by old con- 
straint. 

The partition function of the rigid rotator in the high temperature limit becomes 



Thus, we first evaluate the denominator in the Eq.(12) 

3 3 ^1 

or it is explicitly written as 

Ebf \n q = {1 ^ L Efi - a - ?) w (1 -*>- 
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By defining 

£[1 - (1 - = *, 
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we write it in the following simple form 

3 

We now proceed by using the mathematical relation between (3 and (5' which 
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Substituting Eq. (17) into the Eq. (18), we get 

t' = (2-q) {q+l)/q t l/q . 

From Ref. [9], we have 

By mere substitution, we get 

Z'f\fi) = (2-q) 1/ H- 1/q . 

We also have 



Thus, we formed pf\(3). Next, we use the relation 

Ebf (OT = (V 3) ) 1 - , > (23) 

j 

to obtain the partition function 

= (2 - g)V5 f (l-9)/9. (24) 

Therefore, one may obtain internal energy and heat capacities by using the following relation 
We find 

= (2-q) 1 ' q t 1 / q , (26) 

where /n g a; = gl 1 J" 1 ■ Heat capacity is given by 

C7( 3 ) = ^ r = *(2- (Z ) 1 ^a-«)/«. (27) 

Again, as q — > 1, — ► fc B is obtained. This is the rotational specific heat of the isotropic 
rigid rotator for high temperature limit in accordance with the equipartition theorem. Fi- 
nally, we plot in terms of reduced temperature t=T/9 in Fig. 1. 
ii) Low temperature limit 

At low temperatures, if q is small enough concerning the non-extensive (1-q) part, we 
have T^C 6. We look at first few terms of the summation in the Eq. (3) 

Z f) ^i + 3 [i_2(l-g)^]V(i-9). (28) 

This is nothing but the well-known partition function of the rigid rotator in low temperature 
limit in MB statistics if (1-q) is small enough. This is the case if we make q closer to 1. For 
such a choice, the partition function takes the form 

1 + 3exp(-2#/T). (29) 



The generalized internal energy function for this case becomes 

C/f) = -[1 + 3[1 - 2(1 - g)^] 1 /^)]"^! - 2(1 - q)^-} q l {l - q \ (30) 
a a a 

This is a relatively long expression for the internal energy function but if we look at the 
1 limit again, we see that it is of the form 

u^ i = pi eM _ 2e/n (31) 

with 1/a = ksO in the q— > 1 limit. Calculation of the specific heat reads 

cf = Jljli + 3[i - 2(1 - ,)§] w-.>]-.[i - 2(1 - 

369 : [1 + 3[1 - 2(1 - ,)2]VP-«>]-«-i[i - 2(1 (32) 



a 2 kT 2 a a 

To modify the equations above in accordance with the third constraint, we substitute the 
Eq. (19) into the Eq. (3) and take the first two terms to get the low temperature limit 
partition function 

Zf 1 + 3[1 - 2(1 -q)(2- g )-(9+i)/90V5 T -i/9]i/(i-9). (33) 
Using this partition function, we immediately get the specific heat value as 

Cf = |^(2 - g)-^ +1 ^[l + 3(1 - 2(1 - q)(2 - 5 )-(9+i)/9 f -5)V(i-9)]-5 

[1 - 2(1 -q)(2- q )-^+ l )lH-ifl^-i). (34) 

As q — ► 1, it reduces to 

1 9 

c » s, ^^ exp (- 2 ^' < 35) 



This limiting value of heat capacity is obtained by observing that — ^ — is almost equal to 

Z q^l 

unity for low temperature partition function values. This consideration also holds for Eq. 

(33) 

The plots related to low temperature case are illustrated again in terms of the reduced 
temperature in Fig. 2. The interesting feature in these plots is that specific heat function 
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of the rigid rotator in the low temperature regime attains the same shape as the classical 
one but with a narrower width and a shift in the peak to the left. By increasing q, specific 
heat function attains the same shape as the classical one but with a narrower width and a 
shift in peak to the left. 



III. ISOTROPIC NON-RIGID ROTATOR 



In this case, energy levels are given by [10] 



2// J + lj 2I 2 kR 2 



E^ ri9ld) = £-jU + 1) - T^fU + I) 2 - (36) 



where Iq = /iRq and k is spring constant. We have neglected the third order term j 3 (j + l) 3 
in the above expression. It is better to rewrite equation above in the simple form 

rp(nonrigid) 

= f3 ^ = Bj(j + 1) - Dj 2 (j + I) 2 . (37) 

where D is called centrifugal distortion constant. The first term on the right hand side simply 
corresponds to the rigid rotator part which had been studied in detail in the previous section, 
and non-rigidity term is explicitly given as, apart from a multiplicative factor which is very 
small compared to the rigidity part factor, j 2 (j + l) 2 . The partition function in the NEXT 
formalism reads 

oo 

z (nonrigid) = J-^ + !)[!_ ^ _ (nonrigid)^ 1/(1 _,) _ ( g g) 

j=0 

By substituting the related energy equation to the equation above, we get 

oo 

z (nonrigid) = + 1)[1 - /3(l - q)Bj(j + 1) + (3(1 - q)Dj 2 (j + l) 2 ] 1 /^ . (39) 

j=0 

We again look for its solutions in the high and low temperature limits. 
i)High temperature limit 

We make the same assumptions used in the the rigid rotator case. 

/•oo 

z {non ngl d) = / ^ _ pft _ + ^ _ q ^ Dx 2n/(l-<l)^ ( 4 0) 

JO 



or simply 

Z q = / dxlax 2 + bx + l] 1/{1 - q) . (41) 
Jo 

where 

a = P(l-q)D and b = -(3(l-q)B. (42) 

The integral above can be rewritten in factorial form 

i r°° 1 1 

Z q = {mn)— / dx{x + -) l '^\(x + -WM (43) 
Jo m n 

where m = b — n and n — | [1 ± (1 — 4p-) x / 2 ]. We shall use the following general form to 
evalute the above integral 

[°° dxx v ~\x + PY^x + ^y* = [3-»Y~ (7 Beta(v, fi - v + a^F^y; n + a;l - I), (44) 
Jo p 

where \i = v — 1, 7 = ^ and (3 — ^. 2 Fi(p, /i; 7; x) is the hypergeometric function 
given below 

2 F i(p>nr,x) = T ^ )T ^ } _ dttr-^i-ty-^ii-xt)-", (45) 

where the parameters satisfy z/ > and fj, > v — a. Here we see that z/ = 1 and g > 1. 
Thus, we simply get 

Vl = a = -L, (46) 

By using the condition /i > v — a, we obtain an upper limit for q as q < 3. Therefore, we 
have a solution range for q as 1 < q < 3. Thus we get the partition function as 

Z 9 = jf 1 ^- 2^(1, l/(g - 1); 2/(g - 1); 1 - -). (47) 
(3 — gjn n 

. In order to get the compact expression for the partition function above, we made use of 
the following identity: 

Beta(l, - 1) = q —±. (48) 

q — 1 3 — g 
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The resulting internal energy and heat capacity expressions may be found by using 
Equations (7) and (10) respectively together with the partition function above. Variation of 
heat capacity in the high temperature limit is plotted in Fig. 3. In all the plots related to 
non-rigid rotator, we have used 10.397 and 4.1xl0~ 5 as B and D parameters respectively 
for HC1 molecule [11]. 

ii) Low temperature limit 

As in the case for isotropic rigid rotator in the low temperature limit, we look at first 
few terms of the summation in the Eq. (39) 

Z q = 1 + 3[1 - 25(1 - q)(3 + 4D(1 - q)(3}^. (49) 

Internal energy term resulting from the partition function above is 

U q = Q(B - 2D)[1 + 2(B - 2D)(q - l)0\^[l + 3(1 + 2(B - 2D)(q - l)(3)^]- q . (50) 

The specific heat function of non-rigid rotator for this case is illustrated in Fig. 4, with 
using B and D constants for HC1 molecule, 10.397 and 4.1 xl0~ 5 again. 

IV. RESULTS AND DISCUSSIONS 

We have studied isotropic rigid and non-rigid rotators in the non-extensive statistics 
by extending the previous work on anisotropic rigid rotator [12]. We have calculated heat 
capacity at high and low temperature limits for each rotator. At high temperature, results 
are in well agreement with classical cases. We also observe some interesting features for 
low temperature case. There happens to be shift to the left in the zeros of internal energy 
functions by increasing q as is plotted in Fig. 2. By increasing q from 1 to 3/2, specific heat 
function preserves the same shape as the q=l case but with a narrower width and a shift 
to the left. One of the interesting features of Tsallis statistics is that in the limiting Tsallis 
index, i.e. as q^ 1, one gets the corresponding Maxwell-Boltzmann result at once. This has 
also been true for the partition functions, generalised energy functions and specific heats in 
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the high and low temperature limits of the isotropic rigid rotator. We have concluded that 
study of extensive systems through non-extensive thermostatistics bears importance [13]. 
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FIGURES 

FIG. 1. Specific heat of the rigid rotator as a function of reduced temperature t=T/0 in the 
high temperature limit. 
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FIG. 2. Specific of the rigid rotator heat as a function of reduced temperature t=T/6 in the 
low temperature limit. 




FIG. 3. Specific heat of the non-rigid rotator as a function of temperature in the high temper- 
ature limit. 
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FIG. 4. Specific heat of the non-rigid rotator as a function of temperature in the low tempera- 
ture limit for q= 1.25, 1.50 and 1.75 respectively starting from above. 
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